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Abstract 

. . . In this article, we study the doubly heavy baryon states Hcc, f^cc, Si,;, and nti, in the nuclear 

04 ' matter using the QCD sum rules, and derive three coupled QCD sum rules for the masses, 

vector self-energies and pole residues. The predictions for the mass-shifts in the nuclear matter 
O ! AMh„ = -1.11 GeV, AMn„ = -0.33 GeV, AMh,, = -3.37 GeV and AMn,, = -1.05 GeV 

can be confronted with the experimental data in the future. 
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'Th '' 1 Introduction 

I . In 2002, the SELEX collaboration reported the first observation of a signal for the doubly charmed 

baryon state S+ in the charged decay mode 5+ — > A+if~7r+ in the charm hadro-production 
(-H . experiment (E781) at Fermilab [I], and confirmed later by the same collaboration in the decay 

mode S+ — pD^K~ with the measured mass Ms^^ = (3518.9 ±0.9) MeV [5]. And there have been 
several theoretical approaches to calculate the doubly heavy baryon masses [3] . In Ref. [4] , we study 
^ I the = i"*" doubly heavy baryon states figg and Sqq by subtracting the contributions from 

|0 I the corresponding ^ doubly heavy baryon states with the QCD sum rules, and make reasonable 

predictions for their masses. 
' The QCD sum rules are a powerful theoretical tool in studying the ground state hadrons both 

in the vacuum and in the nuclear matter, and have given many successful descriptions of the 
hadron properties [5]. In Ref.j2], we study the A-type and E-type heavy baryon states in the 
. nuclear matter using the QCD sum rules, and obtain three coupled QCD sum rules for the masses, 

vector self-energies and pole residues in the nuclear matter. In this article, we extend out previous 
works [H |B] to study the properties of the doubly heavy baryon states Sec, ^cc, '^bb and flbb 
in the nuclear matter. The mass-shifts of the heavy mesons, heavy baryons and doubly heavy 
baryons differ greatly from the corresponding ones of the light mesons and light baryons due to 
the appearance of the heavy quarks, the full propagators of the heavy quarks in the nuclear matter 
undergo much slight modifications compared with that of the light quarks, the gluon condensates 
are slightly modified in the nuclear matter. The upcoming FAIR (facility for antiproton and ion 
research) project at GSI (heavy ion research laboratory) provides the opportunity to extend the 
experimental studies of the hadron properties in the nuclear matter into the charm sector [71 , the 
present predictions can be confronted with the experimental data in the future. The ground state 
light-flavor hadrons in the nuclear matter have been studied extensively with the QCD sum rules 
[SlIllIlOl, the heavy quarkonia J ftp, r]c in the nuclear matter have also been studied with the QCD 
sum rules [TT], while the works on the heavy mesons in the nuclear matter are few [T^ . 

The article is arranged as follows: we study the doubly heavy baryon states Sqq and ^qq in 
the nuclear matter with the QCD sum rules in Sec. 2; in Sec. 3, we present the numerical results 
and discussions; and Sec. 4 is reserved for our conclusions. 
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2 The doubly heavy baryons in the nuclear matter with 
QCD sum rules 

We study the doubly heavy baryon states Sqq and ^qq in the nuclear matter with the two-point 
correlation functions n(p), 



Hip) = tJd'xe^P-^^o\T{J,{x)J,iO)}\^o), 

J,{x) = ^^'^Q\{x)C^^Q,{x)^^^-i,qk{x), (1) 

where the i, j, k are color indexes, Q ~ c,b, the C is the charge conjunction matrix, and the |^'o) is 
the nuclear matter ground state. We use the current Jq{x) to interpolate the Sqq, and obtain the 
current Js{x) to interpolate the figg with a simple replacement q ^ s. The correlation functions 
n(p) can be decomposed as 

n(p) = ns{p^,p-u) + iip{p^,p-u)^ + Uu{p^p-u)i, (2) 

according to Lorentz covariance, parity and time reversal invariance [8j[9]. In the limit four- vector 
— (1, 0), the component Ili{p^ ,p ■ u) reduces to Ili{po,p), where i = s,p, u. 
We insert a complete set of intermediate doubly heavy baryon states with the same quantum 
numbers as the current operators Jq{x) into the correlation functions n(p) to obtain the hadronic 
representation ^5 , then isolate the ground state contributions of the doubly heavy baryons Sqq, 
and use the dispersion relation to recast the three components T\i(pQ^p) into the following form: 

n,to,pi . 1 r..^IW-fl, (3) 

27ri J_oo w - po 



where 



_ 

Anp(a.,p) = -2^*^ [5{uj ~ Ep) ~ 5{uj - Ep)] , 

AIljLu,p) = +2m^2R_E2R. [^(^ _ Ep) - <5(c. - Ep)] , (4) 



Ep - + Eh = - E; = ^M^l^ -t-p^, the M|^^, and A|^^ are the 

masses, vector self-energies and pole residues of the doubly heavy baryon states Sqq respectively 
in the nuclear matter. 

We carry out the operator product expansion in the nuclear matter at the large space- like region 
p^ <^ 0, and obtain the spectral densities at the level of quark-gluon degrees of freedom, then take 
the limit = (1,0), and obtain the three components Yii{po,p) of the correlation function n(p) 

urn- 

MP^,P) = J2^niP0,P){On)p., (5) 

n 

where the C^ipo,?) are the Wilson coefficients, (C'n)p„ — (^o|C'n|^o) — (Cn) + pN{On)N at low 
nuclear density in the linear approximation, the (On) and (On) n denote the vacuum condensates 
and the nuclear matter induced condensates, respectively. We can obtain the imaginary parts of 
the QCD spectral densities through the formula 

Ani(w,p) = limite^o [ni(w + ie,p) - ni(w ze,p)] . (6) 
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We match the hadronic spectral densities with the QCD spectral densities, and multiply both 
sides with the weight function (w — Ep)e^T^, perform the integral du, 



dujAIli{uj,p){u - Ep)e-^ , (7) 
to exclude the negative quasi- particle contributions, and obtain the following three QCD sum rules: 

^H««e-S = r dsff^^^{s)e-^, (8) 

Af««M5^^e-S = r dsp^^^{s)e-^, (9) 

■J (2mQ+m,)2 

A|^,«S|«,e-S = r dspl^^{s)e-^, (10) 



3 r»f 
8^ 



da I d(3a(3{l - a - /3)(s - El){bs - 2,EI) 



Hi 



m^Q a,GG 



l-a 



+ -r-T- I da I dp{l - a - p){s - El) 



127r2 



a 



487r2T2 



/ da I dp{l -a- 

'—)p, / da / df3{l - a 



a (3 



1 + 



2T2 



6{s-El) 



13) 



+ 



+ 



rn^asGG_ f-s 
"l67r2^ ^ 

1 ,asGG 



327r2 



/•l — a 

da / dii(\ -a-P) 
da 



1 1 
1 1 



5{s-El) 
5{s-El) 



/3. 



Z + {l+a + p)ml5{s-El) 



Pn 



/ daa{l-a) 5 + {Sm^Q - 2E'q)5{s - E'q) 



37r2 



/■"/ r ^ ^ 

/ daa(l-a) 2 + - 8E'^)6{s - E}^ 



5 ml -El 



2 + - 



T2 



b{s - Eq) + -2 



daa(l — a) 



4 + 



- 12i;2 



3T2 



(5(5 



(11) 
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mgm^ 
' 487r2 



(— — )pAr da 
^ Joe, Jb, 



1-a 



dp 



I PN 



87r2 



/>1 — a 

da / (i/3 

l-a 



a ;9 " 

1 1 

1 1 " 



1 + 



jt2 

5{s-El) 
5{3-El) 



5{s-El) 



167r2 
27r2 



)pN I da I dl3 



m 



1 + ^5{s-E'q) 



daa(l — a) 



3s + 2m^ - 2El^ 



+ 



+ 



{qgs<jGq)p„ 
{qiDoiDoq)p^ 



J 



daa(l — a) 

a, 



5ml 



-%^]5{s-El) 



daa(l — a) 



{qgs(TGq)p^ ["f 



/■«. 

J OLi 



da 



^ mg{qiDoq)pj^ - {qiDQiDQq)p^ /""^ 



67r2 

_;g f ^mg{qU)p« 
^1 7r2 



(5(5 



rp2 



5{s-El) 



5{s-El) 



daa{l — a) 1 + fnQ5{s — Eq) 



(12) 



Of 



{<1U)P. 

2^2 

Pn 



daa(l — a) 



27r2 

{qiDoiDoq)p^ 
7r2 



S + THq 

/ daa(l — a) 



?2 , "^0 



3T2 3T2 



(iaa(l — a) 



8m| 



3r2 



rp2 



5{8-El) 



27r2 



da 



U(s-i?^) 



— £■„ 



mq{qq)pK - ^{qHDQq)^ 
37r2 



/ 



daa{l - a) 2 + TOqJ(s - 



(13) 



where iJg = 



g + ^ +P ,'mQ 

_ l-^l-iml/s _ ami 

a% — o ^ yi — ^„ ™ 



T 3 , -C/| 



Q 



„* _ , ,2 _ „0 



QQ 



a(l-a) ^-f^ ' '"Q ~ a(l-a) ' "/ " 2 ' 

0^, the Shqq ^'^^ continuum threshold 



parameters, the T"^ are the Borel parameters, a^, /3j ^ 0, a/ ^ 1 in the spectral densities where the 
S functions S{s — Eq) and S(s — Eq) appear. We can obtain the masses Mi^^^, vector self-energies 
S| and pole residues in the nuclear matter by solving the three equations Eqs.(8-10) with 
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simultaneous iterations. With the simple replacements of the corresponding parameters, such as 
niq — > nis, {qq)pM ~^ (^'S)pjv7 (9^9) PN ~^ {^''^)pNJ 6tc, we can obtain the corresponding three QCD 
sum rules for the doubly heavy baryon states ^Iqq. 

3 Numerical results and discussions 

The input parameters are taken as {qU)pN = ^Pn, (s"''s)p„ = 0, {qq)pj^ = (qq) + „,°^„,^ Pn, 
{ss)p, = (ss) + yj^^^PN, m„ + - 12MeV, ctn = (45 ± 10) MeV, (^)p„ = {^) - 
(0.65±0.15)GeVpAr, {^) - (0.33 GeV)^, {qUDoq)p„ = O.lSGeVp^v, {snDos)p^ - t^(ss>p« + 
0.02GeVpAr, {qHDoiDoq)p^ + -^{q''gsaGq)p^ = 0.031 GeV^JV, {sHDoiDos)p„ + -^{s'' gscrGs) p^ = 
yO.031 GeV^PAT, (qgsCrGq) ml{qq), (sgsCrGs) = mg(ss), {qiDoiDoq)p„ + l{qgs<jGq)pj^ = 0.3GeV^pjv, 
(sii:)o«£'os)p„ + |(s.gsCrGs)p„ = yO.SGeV^pAr, {qgsaGq)p^ = {qgs(TGq)+3.0GeV^ pN, {sgsaGs)p^ = 
(sgscrGs) +y3.0GeV'^PN, {q^gs^Gq)p„ = -0.33 GeV^pjv, {s^gs(7Gs)p^ = -yO.33 GeV^PAr, (gg) = 
-(0.23GeV)3, (ss) = 0.8(qq), = 0.8 GeV^ pat = (O.llGeV)^, y = 0.3 ± 0.3, m, = 6MeV, 
rus = 0.13 GeV, rUc = 1.35 GeV and mi, = 4.7 GeV at the energy scale fi = 1 GeV 

We can recover the QCD sum rules in the vacuum by taking the limit p^v = 0, then differentiate 
the Eqs.(8-9) with respect to yj- respectively, and eliminate the pole residues X^qq (here we smear 
the asterisk * to denote the pole residues in the vacuum), and obtain two QCD sum rules for the 
masses M^^^ with respect to the spinor structures ^ and 1, respectively. On the other hand, we 
can divide Eq.(9) by Eq.(8) to obtain the QCD sum rules for the masses AI^qq- From Eq.(lO), 
we can see that if we take the limit pjv = 0, "^^qq 7^ 0, due to appearance of the term mq{qq) p^^, 
however, such terms arc small enough to be neglected safely. The QCD spectral densities from the 
operator product expansion in the vacuum and in nuclear matter by taking the limit pN = have 
slight discrepancies, the discrepancies cannot result in remarkable differences on the masses and 
pole residues. 

In Ref . [Ij , we study the doubly heavy baryon states Sqq and ^qq by subtracting the possible 
contributions from the corresponding negative-parity doubly heavy baryon states with the QCD 
sum rules. In that reference, we choose the spinor structure 1+7", take the Borel parameters as 
^ (2.8-3.8)GeV2, T^^^ = (3.0-4.0) GeV^, Ti^^ = (7.7-9.1) GeV^, T^^^ ^ (7.9-9.3) GeV^ 
the continuum threshold parameters as s|^^ = (4.2 ± 0.1 GeV)^, s^^^ = (4.3 ± 0.1 GeV)^, s| = 
(10.8 ± 0.1GeV)2, s^^^ = (10.9 ± O.lGeVT^, and obtain the masses Ms,, = (3.57 ± 0.14) GeV, 
Mn,^ = (3.71±0.14)GeV, Ms,, = (10.17±0.14) GeV, Mo,, = (10.32±0.14) GeV. The predictions 
are consistent with the value from the SELEX coUaboration M^^^ = (3518.9 ±0.9) MeV [Tl[2]. For 
the technical details, one can consult Ref.j4j. 

In this article, we choose the spin structures 1,]^ and ji, see Eq.(2), take the continuum threshold 
parameters as sl^^ = (17.5 ± 0.5) GeV^ s^_^^ = (18.0 ± 0.5) GeV^ s|^^ = (117.0 ± 1.0) GeV^ and 
^%.bb ~ (119-0 i 1.0) GeV^ consulting Ref.|4j, and vary the Borel parameters to reproduce 
almost the same masses as Ref. [4]. In calculations, the Borel parameters are taken as T~ = 
(2.4 - 2.8) GeV^ T^^^ = (2.4 - 2.8) GeV^ r|^^ = (7.0 - 7.5) GeV^ and T^^^^ = (6.7 - 7.2) GeV^ 
respectively. 

The interpolating currents Jq {x) have nonvanishing couplings to both the positive- and negative- 
parity doubly heavy baryon states. For the QCD sum rules in the nuclear matter, we cannot choose 
the spinor structure 1 -I- 7° to avoid the possible contaminations from the negative-parity doubly 
heavy baryon states. In this article, we choose the spinor structures 1 and p in stead of the 
spinor structure 1 -t- 7*^ in the limit pjv = 0, the ratios between the negative- and positive-parity 
contributions are about —(0 — 9)% and —(4 — 8)% in the charm and bottom sectors, respectively 
[4]. For the technical details, one can consult Ref. [4]. The contaminations from the negative-parity 
doubly heavy baryon states are rather small. We can also take into account the contributions 
from the negative-parity doubly heavy baryon states explicitly as Ref.[13j. and perform detailed 
analysis, the tedious analysis maybe our next work. 
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If we take the limit triu = ma = mq — >■ in the QCD sum rules for the nucleons, the hadronic 
parameters M'^ and S]^ in the nuclear matter can be approximated as 

S]V = ^(g+g)p., (14) 

respectively. Such relations are simple and elegant, however, the predictions change monotonously 
with variations of the Borel parameter T^. In the present case, if we take the limit ruq — )• and 
neglect other terms of minor importance, we can obtain the following relations, 

= g{T^mQ,sg^^){qU)p., (15) 

from Eqs.(8-13). The formal notations f (T-^ , mn , s'i ) and g{T'^,mQ,s% ) are complex func- 
tions of the variables T^, mq and sSgg- It is difficult to carry out the integrals over the variables 
s, a, j3 analytically due to the appearance of the heavy quark mass mg, we expect that the 
f{T^, mn, s~ ) and giT"^, mn, s~ ) arc sensitive to the Borel parameter T^, the numerical cal- 
culations confirm such conjecture. 

Finally, wc obtain the masses and pole residues in the vacuum Mh^^ = (3.53 ± 0.32) GeV, 
Mn^^ = (3.61 ±0.29) GeV, Ms,, = (10.14±0.47) GeV, Mn,, = (10.27±0.46) GeV, Ah„ = (0.117± 
0.046) GeV, Xn,^ = (0.128 ± 0.047) GeV, Ah,, = (0.287 ± 0.162) GeV, A^,, = (0.359 ± 0.211) GeV; 
the masses, pole residues and vector self-energies in the nuclear matter AI~ = (2.43 ± 0.27) GcV, 
= (3.28 ± 0.26) GcV, M£^^ = (6.77 ± 0.38) GcV, Mj*j^^ = (9.22 ± 0.41) GcV, Ai = (0.033 ± 
0.005) GeV, A* = (0.078 ± oTo21) GcV, Al^^ = (0.008 ± o'.()02) GcV, A*j^^ = (0.075To.033) GeV, 
= (0.139 ± 0.005) GcV, = -(0.olj6 ± 0.001) GcV, = (0.507 ± 0.019) GoV, E^j^^ = 

(076b9 ± 0.001) GcV. The mass-shifts in the nuclear matter are AMs,, = -(1.11 ± 0.04) GeV, 
AMn,, = -(0.33 ±0.03) GeV, AMh„ = -(3.37 ±0.09) GeV, AMq,, = -(1.05 ±0.06) GeV, where 
the mass-shifts AM defined by AM = M* — M are the scalar self-energies. The uncertainties 
come from the Borel parameters T^. 

The mass modifications in the nuclear matter are about —31%, —33%, —9% and —10% for the 
doubly heavy baryon states Sec, '^bb, ^cc and flbb, respectively, which are consistent with the quark 
condensates modificatfons in the nuclear matter ^^"^^"^^"^^"^ « -41%, ^"'^"^^"^^"'^ « -(15 ± 15)%, 
where the uncertainty ±15% comes from the uncertainty of the parameter 6y = ±0.3. There appear 
additional quark condensates associated with the light flavor quarks, such as {q^q)pj^, {q\Doq)p^ , 
(qHDoiDoq) , (qiDoiDoq) , etc, their contributions are not large but cannot be neglected safely. 
The Sqq and ilgg baryons have two heavy quarks besides a light quark, the heavy quark interacts 
with the nuclear matter through the exchange of the intermediate gluons, and the modifications 
of the gluon condensates in the nuclear matter are mild, {^^^)pr, « 0.93(^i^) , we expect that 
the mass modifications are smaller than that of the nucleons, which have three light quarks, uud 
or udd, and the approximation in Eq.(15) makes sense. 

If we take into account the uncertainties of the continuum threshold parameters, Ss~ = 
±0.5GeV^ Ss'i^^^ = ±0.5GeV^ Ssg^^ = ±1.0GeV^ ,5s^^^ = ±1.0GeV^ additional uncertainties 
(5AMh„ = ±0.03 GcV, SAMn,^ ='±0.01 GeV, SAM^.l = ±0.10 GeV, SAMn,, = ±0.04 GeV 
for the mass-shifts AM are introduced. The uncertainty 5y = ±0.3 leads to the uncertainties 
5AMn,, = ±0.43 GeV and (JAMq,, = ±1.13 GeV. The uncertainty SaN = lOMeV leads to the 
uncertainties (5AMh„ = ±0.37 GeV, SAMq^^ = ±0.10 GeV, ^AMh,, = ±1.05 GeV and ^AMn,, = 
±0.27 GeV, respectively. 

We can also take into account the perturbative 0{a^) corrections in the leading logarithmic 
approximation through multiplying the QCD spectral densities by the anomalous-dimension fac- 
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tors, 



(16) 



where the anomalous-dimensions of the interpolating currents are T j — , and the anomalous- 



dimensions of the local operators O are r^^^^ — -gj-rfjvj' ^"ig{??> ~ ^(s?gg> = 0, and the Nf is the 
flavor number. The anomalous-dimension factors can lead to the uncertainties — (13— 17)%, 

^ = (10 - 13)%, ^ - (28 - 36)%, ^ = (21 - 27)%, ^ (16 -21)%, ^ = 

(10 - 13)%, ^-^P^ = (34 - 46)%, = (20 - 26)% for the values Kqcd = (200 - 300) MeV. 

The masses in the vacuum and in the nuclear matter are increased significantly, we can normalize 
the masses in the vacuum to the values M^^^ = 3.53 GeV, Mq^^ — 3.61 GeV, Mh^^ — 10.14 GeV, 
Mubb = 10.27 GeV, and estimate the uncertainties of the mass-shifts due to the energy scales as 
^ = -(5 - 7)%, '4mt = (4 - 5)%, ^ = -(10 - 14)%, ^ = (6 - 8)%. 

We can refit the Borel parameters and threshold parameters to reproduce the experimental 
data approximately the Borel parameters (or the effectively energy scales) and the threshold 
parameters have some correlations. In calculations, we observe that larger threshold parameters 
can lead to smaller masses and cancel out the enhanced factors induced by the non-zero anomalous- 
dimensions in Eq.(16), however, the contributions from the high resonances and continuum states 
are included in. So we can fix the threshold parameters and take larger Borel parameters to cancel 
out those enhanced factors, see Eqs. (14-15), then the uncertainties of the mass-shifts induced by 
the energy scales are about a few percents, and can be neglected safely. On the other hand, we can 
introduce Borel parameter dependent threshold parameters and study the systematic uncertainties 
[l4] , it is a hard work before the experimental data are enough and the precise values of the pole 
residues are known, we postpone those works in the future. 

If we take the loffe current to interpolate the proton, the QCD sum rules indicate that there 
exists a positive vector self-energy S]^ = (0.23 — 0.35) GeV with the typical values of the relevant 
condensates and other input parameters and a reasonable negative scalar self-energy with the 
suitable parameters There exists substantial cancelation between the scalar and vector self- 
energies, the self-energies and T,"^, which correspond to the real energy- independent optical 
potentials S and V, satisfy the relation E|^f/E}(f « —1 in the leading order approximation. While 
the mean-field models predicate that the typical self-energies of the nucleons in nuclear matter 
saturation density are « —350 MeV and E](, w -|-300MeV respectively, the effective non- 
relativistic central potentials S + V are about tens of MeV. In the present case, E2. -|- EH, = 
-0.97 GeV, E^ E^ =-0.34 GeV, El + EH = -2.86 GeV, E?, -h EH, = - L04 GeV,"the 
net optical potentials S + V are large, as the self-energies |E^| ^ jE'*!. The present prediction of 
the mass-shift SAI^^^ = —1.11 GeV can be confronted with the experimental data from the CBM 
and PANDA collaborations in the future [7] , where the properties of the charmed baryons in the 
nuclear matter will be studied. The Sec have interesting properties, for example, the theoretical 
predictions of their lifetimes based on different quark models are in agreement with each other 
|15j . but about one-order larger than the upper limit of the experimental data [1]. The precise 
measurement of the lifetime by the PANDA and LHCb collaborations in the future maybe shed 
light on the apparent discrepancy. 



4 Conclusion 

In this article, we extend our previous works on the A-type and E-type heavy baryon states to study 
the doubly heavy baryon states Sqq and i^QQ in the nuclear matter using the QCD sum rules, and 
derive three coupled QCD sum rules for the masses, vector self-energies and pole residues in the 
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Figure 1: The masses from the QCD sum rules in the vacuum and in the nuclear matter versus the 
Borel parameter T^, the A, B, C and D denote the Sec, ^cc, '^bb and flbb baryons, respectively. 
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nuclear matter, then take the limit pjv = to recover the QCD sum rules in the vacuum, finally 
obtain the values of the masses and pole residues in the vacuum, and the masses, vector self-energies 
and pole residues in the nuclear matter. The numerical results indicate that the mass-shifts in the 
nuclear matter are about AMh,, = -1.11 GeV, AAfo^^ = -0.33 GeV, AMh,, = -3.37 GeV and 
AMn^b = —1.05 GeV, respectively, which can be confronted with the experimental data in the 
future. 
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